Theoretical 2D Raman band of strained graphene 



Faculty of Physics, 



Valentin N. Popov 

University of Sofia, BG-1164 Sofia, Bulgaria 



o 



Philippe Lambin 

Research Center in Physics of Matter and Radiations, 
University of Namur (FUNDP), B-5000 Namur, Belgium 
(Dated: March 27, 2013) 

We study the 2D Raman band of in-plane uniaxially strained graphene within a non-orthogonal 
tight-binding model. At non-zero strain, the obtained 2D band splits into two subbands at strain 
angles 0° and 30° or into three subbands at intermediate angles. The evolution of the 2D subbands 
is calculated systematically in the range of the accessible strains from —1% to 3% and for the 
commonly used laser photon energy from 1.5 eV to 3.0 eV. The strain rate and dispersion rate 
of the 2D subbands are derived and tabulated. In particular, these two quantities show large 
variations up to 50%. The results on the 2D subbands can be used for detecting and monitoring 
strain in graphene for nanoelectronics applications. 
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I. INTRODUCTION 

The successful isolation of flakes of the two- 
dimensional crystal graphene and the measured high 
thermal and carrier mobility open new possibilities for 
application in nanoelectronics 1 and nanophotonics^ It is 
important to investigate the state of strain in graphene, 
because strain can effectively modify the electronic struc- 
ture and phonon dispersion,- and thus change the trans- 
port properties of this material A This can be done eas- 
ily and in a non-destructive way by means of Raman 
spectroscopy^ 

The Raman spectrum of perfect graphene is dominated 
by an intense first-order Raman band, the so-called G 
band, and a few second-order Raman bands£ The most 
intense one among the latter is observed at w 2700 cm -1 
and is termed the 2D band. This band is dispersive, i.e., 
it depends on the laser photon energy El- It is excited 
by double-resonant scattering processes^ favored by the 
specific cone-like electronic structure near the Fermi en- 
ergy. Raman measurements on strained graphene have 
been performed at different strain magnitude, strain di- 
rection, and laser photon energy&~— In all these studies, 
only one or two peaks of the 2D band could be resolved 
and their shift rates have been derived. In two similar 
measurements) 10 ' 12 the determined strains differ by a fac- 
tor of 2, which shows that the calibration of strain can 
still be an issue. 

The dependence of the shift rate on strain can be ob- 
tained by calculation of the 2D band under various strain 
conditions. The complexity of these calculations at the 
ab-initio level hinders the detailed investigation of the 
dependence of this band on strain and laser photon en- 
ergy. The theoretical work has so far been limited to the 
calculation of the 2D peak position^ - — while the in- 
tensity of the peaks is either treated approximately; 1 b 12 
or, in a limited number of cases, is calculated from first 
principles.— Despite the success in modeling the 2D Ra- 
man band of graphene, a detailed investigation of its be- 
havior under strain is mandatory for monitoring strain 



in this material by means of Raman spectroscopy. 

Recently, we have reported the predicted 2D band of 
graphene at strain magnitude e = 1% and El = 2.5 eV 
within a non-orthogonal tight-binding (NTB) models 
The very good agreement of our results with the avail- 
able experimental data supports the applicability of this 
model to the 2D band of strained graphene. 

Here, we use the NTB model for large-scale calcula- 
tions of the 2D band in graphene under uniaxial strain 
in the range of the accessible strain magnitudes and dis- 
tinct strain directions, as well as for the common laser 
photon energies. The NTB model is presented in Sec. II. 
The obtained results are discussed in comparison with 
available experimental data in Sec. III. The paper ends 
up with conclusions (Sec. IV). 



II. COMPUTATIONAL DETAILS 

The NTB models uses Hamiltonian and overlap ma- 
trix elements obtained from an ab-initio study^ This 
model has no adjustable parameters, which is an advan- 
tage over the majority of tight-binding models. It also 
allows one to determine the total energy and forces on 
atoms. This feature is necessary for the relaxation of 
the atomic structure of graphene and for the calculation 
of the phonon dispersion. The latter is performed us- 
ing a perturbative approach within the NTB model with 
electron-phonon matrix elements derived within the same 
modelJ^ The structure relaxation and phonon disper- 
sion calculations require the summation over the Bril- 
louin zone of graphene. The phonon frequencies con- 
verge to within 1 cm -1 on increasing the size of the 
Monkhorst-Pack mesh of points up to 40 x 40. The calcu- 
lated high-frequency phonon branches are overestimated 
by the NTB models However, very good agreement with 
experiment can be reached by scaling the phonon fre- 
quencies of these branches by a factor of 0.9j^ For this 
reason, the phonon frequencies are systematically scaled 
by a factor of 0.9 for all calculations of the 2D band in 
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FIG. 1. (a) The hexagonal atomic structure of graphene. The 
atoms are depicted by solid circles and the solid lines are the 
interatomic bonds. The strain direction is defined by the an- 
gle 6 relative to the x axis, which is chosen along a zigzag 
of carbon bonds. The limiting cases of Z-strain and A-strain 
are also shown, (b) Schematic representation of a double- 
resonant Raman scattering process, which contributes to the 
2D band. An electron is scattered between the cone-like con- 
duction bands at the K and K points of the hexagonal Bril- 
louin zone of graphene by a K point phonon with a wavevec- 
tor, denoted by an arrow. The electron is scattered back by 
a phonon with an opposite wavevector. The three scattering 
paths KK are denoted by 1, 2, and 3. Similarly, holes can be 
scattered by phonons between cone-like valence bands at the 
K and K points. 



the paper. 

The Raman intensity of the 2D band is calculated us- 
ing fourth-order terms in quantum-mechanical perturba- 
tion theory.— The required electron-photon and electron- 
phonon matrix elements, and electronic linewidth are ob- 
tained within the NTB model. The expression for the in- 
tensity contains summation over all electronic states and 
over all phonons in the entire Brillouin zone. The conver- 
gence of the Raman shift of the 2D band within 1 cm -1 
is achieved by increasing the size of the Monkhorst-Pack 
mesh of points up to 800 x 800. 



III. RESULTS AND DISCUSSION 

We calculate the electronic structure and phonon dis- 
persion of the 2D band of uniaxially strained graphene 
at strain magnitude e = —1%, 1%, 2%, and 3% and 
strain direction, defined by the angle 9 relative to a 
zigzag of carbon bonds, with values 9 = 0°, 10°, 20°, 
and 30° (Fig. la). The interval from 9 = 0° (Z-strain) 
to 9 = 30° (A-strain) covers all distinct strain directions 
in graphene. The electronic and phonon dispersions of 



strained graphene have been discussed a lot in the litera- 
ture (see, e.g., Ref.ll). Previously, we have demonstrated 
that both dispersions are correctly described by the NTB 
models and for this reason, here, we focus on the study 
of the 2D band and its dependence on strain. We con- 
sider laser photon energies El =1.5, 2.0, 2.5, and 3.0 eV, 
and restrict ourselves to parallel and perpendicular light 
polarizations relative to the strain direction for both po- 
larizer and analyzer. Finally, we are concerned only with 
the Stokes part of the Raman spectra. 

In order to facilitate the understanding of the following 
results, we describe briefly the double-resonant processes. 
In such a process, an electron (hole) is scattered be- 
tween states of the cone-like conduction (valence) bands 
at the K and K points of the Brillouin zone of graphene 
(so-called Dirac cones) by two phonons with opposite 
wavevectors close to the K point (Fig. lb) . The resulting 
2D band is centered at roughly twice the frequency of the 
phonons. Each K point has three neighboring K points 
and therefore there are three scattering paths KK . Due 
to the hexagonal symmetry of perfect graphene, the scat- 
tering along the three paths yields equal contribution to 
the Raman shift and intensity of the 2D band. 

It has been debated in the literature about the rela- 
tive contribution to the 2D band from phonons close to 
the TK and KM directions in the Brillouin zone»i£ Such 
phonons and the associated processes are often called 
"inner" and "outer", respectively. Some authors have 
concluded that the contributions from these processes 
are comparable ; 11 ^ while others have found that the 
contribution from "outer" processes is negligible com- 
pared to that from "inner" ones.— ^ It has also been ar- 
gued that the phonon wavevector for a seemingly "outer" 
process can be transformed by a reciprocal lattice vec- 
tor to a phonon wavevector for an equivalent "inner" 
process, "thereby dismissing the notion of inner and 
outer processes" In a recent study, we have proposed 
a quantitative criterion for distinction between the two 
processes.— We have noted that (1) the phonons with 
largest contribution to the 2D band belong exclusively 
to the TO phonon branch and (2) the TO branch is the 
fifth one close to the TK direction and the sixth one close 
to the KM direction (in order of increasing frequency). 
Based on these observations, we have defined unambigu- 
ously "inner" and "outer" processes as processes due to 
phonons from the fifth and sixth phonon branches close 
to the K point, respectively. Using this definition, we 
have shown that in unstrained graphene the ratio of the 
contributions to the 2D band from "inner" and "outer" 
processes is about 10. 

The calculated 2D band of unstrained graphene has 
a symmetric shape (Fig. 2, dotted line). This result is 
a direct consequence of the dominant contribution from 
"inner" processes and is in agreement with the observed 
symmetric shape of the 2D band^ 
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FIG. 2. Calculated 2D band at = 0° (Z-strain), 10°, 
20°, and 30° (A-strain), at E L = 2.0 eV, and for paral- 
lel/perpendicular light polarization (solid/dashed lines) in 
comparison with that at zero strain (dotted line). In most 
cases, two Raman features are clearly seen. They correspond 
to the experimentally observed 2D~ and 2D + bands. 



A. Strain dependence of the Raman intensity of 
the 2D subbands 

The calculated 2D Raman spectra at different strain 
magnitude and at El = 2.0 eV, for parallel and perpen- 
dicular light polarization are shown in Fig. 2. The spectra 
at El = 1.5, 2.5, and 3.0 eV are also calculated but are 
not provided here. In all cases, two or three 2D subbands 
are present. In order to explain their origin, we note that 
in strained graphene, the electronic and phonon disper- 
sions are modified differently along the three paths KK . 
Therefore, the 2D band splits into two or three subbands 
depending on the strain direction. The subbands, com- 
ing from paths 1, 2, and 3, are denoted here by 2Di, 2D2, 
and 2D3 , respectively. These subbands are mainly due to 
"inner" processes, while the contributions from "outer" 
processes give rise to minor peaks and kinks, which de- 
form slightly the shape of the 2D subbands. 

The calculated spectra for Z-strain graphene are given 
in Fig. 2 (upper left panel). In this case, scattering paths 
2 and 3 are equivalent (see Fig. lb). Under tensile strain 
(e > 0) the 2D subbands are red-shifted, while under 
compressive strain (e < 0) they are blue-shifted. The red 
(blue) shift increases with increasing strain (in absolute 
value). There are only two distinct 2D subbands: the 
composite 2D3+2D2 band and the 2Di subband. The 
splitting between them increases with increasing e. For 
parallel polarization, only the composite band is clearly 



seen, while, for perpendicular polarization, both sub- 
bands can be observed in the Raman spectrum. For ten- 
sile strain, the composite band has smaller Raman shift 
than the 2Di subband and vice versa for compressive 
strain. 

For the intermediate values of the strain direction an- 
gle 8 = 10° and 20°, there are three 2D subbands, as seen 
in Fig. 2 (upper right and lower left panels). The separa- 
tion between the 2D3 and 2Di subbands increases with 
increasing e and decreases with increasing 8. At 8 = 10° 
and e > 0, and 8 = 20° and e < 0, the 2D 3 and 2D 2 sub- 
bands overlap considerably and under such strains and 
angles only two 2D peaks can be observed: one peak, 
corresponding to the overlapping subbands 2D3 and 2D2, 
and another peak, corresponding to the subband 2Di. 
Figure 2 (lower right panel) shows that for 8 = 30° only 
two Raman features are present, namely, the 2D3 sub- 
band and the composite 2D2+2D1 band. The red (blue) 
shift of the subbands changes with e similarly to the cases 
of the values of 8 considered above. 

Next, we compare the shape of the obtained 2D band 
to available experimental Raman spectra of strained 
graphene. As a general feature, two 2D Raman peaks are 
usually observed. They are termed 2D - and 2D + band, 
the latter having larger shift than the former J£r— Mea- 
surements on Z- and A-strain graphene at tensile strains 
up to 3.2% and 3.8%, respectively, have been performed 
at El — 2.33 eV (Ref. VWt) . The measured Raman spec- 
tra show mainly one intense peak and a minor peak or 
a kink, especially, for Z-strain at large e. Such a behav- 
ior is qualitatively reproduced by our spectra in Fig. 2 
(upper left and lower right panels). The asymmetry of 
the 2D band shape, namely, higher minor peak for per- 
pendicular polarization than for parallel polarization, is 
also predicted here. In another experiment Raman 
measurement have been performed on a graphene flake 
at 6 = 18.4° and tensile strain e up to « 0.62 (0.58)%, 
and on another flake at 6 = 10.8° and e up to w 1.17 
(1.17)% for parallel (perpendicular) light polarization; in 
both cases El = 1-58 eV. In neither case, the 2D band 
could be fitted with more than two Lorentzian functions, 
because of the small 2D2— 2D 3 band splitting at such rel- 
atively small El- Our spectra in Fig. 2 are in accord with 
these observations. In Ref. |Tl|, the overall symmetric 2D 
band of unstrained graphene has been fitted with two 
Lorentzian functions, attributed to "inner" and "outer" 
processes. It has been concluded that the two types of 
processes give comparable contributions to the band. In 
a recent paper ^ close to Z- and A-strain graphene sam- 
ples have been studied by Raman spectroscopy at strains 
up to 1.19 (1.15)%, at E L = 2.41 eV and for parallel 
(perpendicular) light polarization. By analyzing the 2D 
band shape and splitting, it has been argued that this 
band comes mainly from "inner" processes. Our simu- 
lated spectra exhibit similar 2D band splitting and shape, 
as the latter experimental ones, in support of the rela- 
tively small contribution from "outer" processes. 

The 2D Raman bands can be characterized by their 



4 



TABLE I. Calculated shift rates of the subbands 2Di and 2D3 (in cm~ 1 /%), corresponding to the experimentally observed 
subbands 2D + and 2D - , in comparison with available experimental and ab-initio data. The relative deviation of the NTB 
values from the experimental (ab-initio) ones is given in the columns labeled "Rel. dev." The second row for each couple of 
values for E_l and 9 contains the NTB values scaled by 0.8 and the corresponding relative deviation (see text). 
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linewidth and Raman shift. The calculated full width at 
half maximum increases steeply with El (in eV) at e = 
as 13.2 + 0.067££ (in cm" 1 ). The predicted linewidth 
of 17.8 cm -1 is more than twice smaller than the mea- 
sured one of as 40 cm" 1 at E L = 2.33 eV (Ref.[l3). Our 
linewidth of 15 cm" 1 is by 50% smaller than the averaged 
one of 24 cm -1 for two samples with 9 w 20%, measured 
at E L = 1.58 eV (Ref.dD). At E L = 2.41 eV, we calculate 
the value of 18.6 cm" 1 , which is about 20% smaller than 
the experimental one of ~ 22 cm" 1 (Ref. 12). The dis- 
agreement of the theoretical linewidth with the observed 
one can be due to electron scattering mechanisms, un- 
accounted for here, or to any extrinsic mechanism that 
reduces the lifetime of the phonon. 

The theoretical linewidth increases weakly with e at 
the rate of a few cxcT 1 /% at E L = 2.0 eV (see Fi& 2). 
This behavior corresponds to that reported in Ref. [l0| and 
H3 . In Ref. [Tl|, the measured 2D band at tensile strain 
has been fitted by a single Lorentzian and the obtained 
linewidth has been found to increase steeply with e. In 
this case, the linewidth is roughly equal to the sum of the 
halfwidths of the lower and upper 2D subbands, and the 
separation between them. Our estimated separation of 
11 cm" 1 corresponds well to the increase of the measured 
linewidth of ~ 12 cm -1 at e 0.5%. 



B. Dependence of the Raman shift of the 2D 
subbands on e 

First, we compare the calculated red shift of the 2D 
band to the available experimental data. One of the 
early measurement of these shifts has been performed 
on graphene flakes in polymer films using two- and four- 
point bending setups at E L = 2.41 eV (Ref. @). No split- 
ting of the 2D band has been observed and the red shift 
per unit strain (so-called shift rate) of —64 cm~ 1 /% has 



been obtained. A three-point setup for Raman measure- 
ments on polymer-imbedded graphene flakes has been 
used for systematic investigation of A- and Z-strain sam- 
ples at E L = 2.33 eV. 10 The shift rates of the two 
observed peaks of the 2D band are —16.3 and —29.7 
cm" 1 /% for Z-strain and —21.7 and —30.5 cm" 1 /% for 
A-strain. These values are at least twice smaller than 
those of Ref. |9j. Later on, the analysis of the Raman 
data of strained graphene at 9 = 18.4°, measured at 
El = 1-58 eV, has allowed to derive the value of w —30 
cm" 1 /% (averaged) (—46.8 cm _1 /%) for the shift rate 
of the higher (lower) 2D component The shift rates of 
—49.4 and —57.4 cm" 1 /%, calculated from first princi- 
ples, overestimate the experimental values by more than 
20% (in absolute value). In another paper, the pairs of 
shift rates —44.1 and —63.1 cm" 1 /% (near- A-strain) and 
—26.0 and —67.8 cm" 1 /% (near-Z-strain), measured at 
El = 2.41 eV, have been reported^ The ab-initio cal- 
culations for the two samples have yielded for the pairs 
of shift rates —43 and —70 cm~ 1 /%, and —24 and —66 
cm" 1 /%, respectively. These values reproduce well the 
shift rates of the 2D subbands, except for the lower sub- 
band at A-strain, for which the predicted rate is larger 
(in absolute value) than the measured one by about 10%. 

Table I shows that the NTB shift rates are system- 
atically larger than the available experimental data, the 
overestimation ranging from 19% to 173%. Similarly to 
the case of the G band of strained graphene^ the larger 
shift rates of the 2D subbands can be reduced and agree- 
ment with experiment can be reached by scaling them 
by a single factor. For the determination of the scaling 
factor it is reasonable to choose experimental data with 
observed 2D band splitting and available first-principles 
predictions. Moreover, the reliability of the scaling factor 
will depend on the agreement between the experimental 
and theoretical data. Following these arguments, we use 
the experimental shift rates of Ref. [12J because they agree 
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FIG. 3. Scaled Raman shift of subbands 2Di and 2D3 (sym- 
bols) as a function of the strain e at 6 = 0° and 30°, and at 
different values of El- 



well with the ab-initio ones. From the average overesti- 
mation of the data from the latter paper, we derive a 
scaling factor of ps 0.8. The value of this factor can be 
obtained more accurately by using more precise exper- 
imental and theoretical data. In this paper, we accept 
this scaling factor, and the red (blue) shift will be scaled 
by 0.8 everywhere below, in addition to the scaling of 
the phonon frequencies by 0.9. For brevity, the Raman 
shift with scaled red (blue) shift will be referred to as 
the scaled Raman shift. As evidenced from Table I, our 
scaled shift rates overestimate the experimental data of 
Ref. [13 and Ref. [0 in average by « 100% and « 30%, 
respectively. 

The scaled Raman shifts of subbands 2Di and 2D3, 
derived from our spectra, are presented in Fig. 3 as a 
function of e at different values of E L and 8. In all cases, 
except for the 2Di subband at 9 = 0°, the shifts are 
quasi-linear in e. The shift rate is approximated here 
by the slope of the linear curve fitted to the data points 
in Fig. 3. We note that, in all cases, a linear curve is 
sufficient to fit the data except for the 2Dx subband at 
6 = 0°, where a polynomial of at least second degree 
is needed for a satisfactory fit. The shift rates, derived 
from Fig. 3, are presented in Table II. The shift rate 
for intermediate values of El and 9 can be obtained by 
interpolation. 

Finally, we calculate the Raman spectra of Z- and A- 
strain graphene at El = 2.33 eV and light polarization 
as for the measured spectra in Fig. 3 of Ref. [Hi Namely, 
the incident laser light is polarized parallelly or perpen- 
dicularly to the strain direction and the scattered light 
is observed at an angle of 0°, 45°, or 90° relative to the 



TABLE II. Shift rates (in cm _1 /%) of the subbands 2Di and 
2D3 at different values of El (in eV) and 8, derived from 
Fig. 3 
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FIG. 4. Calculated Raman spectra at Z-strain (e = 1.9%) and 
A-strain (e = 1.3%), at El = 2.33 eV and parallel or perpen- 
dicular incident light polarization. The scattered light is po- 
larized at an angle of 0°, 45°, or 90° relative to the incident 
light polarization. The obtained spectra are in qualitative 
agreement with those reported in Fig. 3 of Ref. 03 



incident light polarization. Since our strain rates over- 
estimate the experimental ones of Ref. [13 by « 100%, 
our calculations are performed at twice smaller strain 
magnitudes than the reported ones, i.e., at e = 1.9% (Z- 
strain) and e = 1.3% (A-strain). The obtained spectra 
(Fig. 4) agree qualitatively with the experimental ones. 
Major discrepancy is observed at Z-strain, where the pre- 
dicted 2Di-2D 3 splitting clearly resolved in perpendicu- 
lar geometry is ~ 50% larger than the experimental one. 
Obviously, more experimental and theoretical efforts are 
necessary for the unambiguous determination of the shift 
rates of the 2D band. 
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TABLE III. Dispersion rate (in cm 1 /eV) of subbands 2Di 
and 2D3 at different e (in %) and 9, derived from Fig. 5. 
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C. Dependence of the Raman shift of the 2D 
subbands on El 

The scaled Raman shifts of subbands 2Di and 2D3, 
derived from our spectra, are shown in Fig. 5 as a function 
of El at different e and 9. In all cases, except for the 
2D 3 subband at 9 = 30°, this function is quasi-linear and 
can be characterized by its first derivative, i.e., the slope 
of the corresponding curve. This quantity describes the 
dispersion of the 2D band and is sometimes referred to 
as the dispersion rate. It is approximated here by the 
slope of the linear curve fitted to the data. We note that 
in the case of subband 2D3 at 9 — 30°, a second-degree 
polynomial fits the data more accurately. 

Table III presents the dispersion rate of the 2D sub- 
bands, derived from Fig. 5, in comparison with that for 
unstrained graphene of 89 cm /eV (Ref. 14). In the in- 
terval of the studied tensile strains, the dispersion rate 
of the 2Di subband increases up to 20%. The dispersion 
rate for the 2D3 subband at 9 = 0° and 10° increases 
slightly, whereas it shows a decrease down to 20% at 
9 = 20° and by almost 50% at 9 = 30°. Therefore, the 
variation of the Raman shift with E L is not negligible and 
has to be taken into account for the correct assignment 
of the 2D components of strained graphene. 

The dependence of the dispersion rate of the 2D band 
of strained graphene on E L has been investigated experi- 
mentally in the case of Z-straini^il The reported rates are 
93, 99, and 106 cm~ 1 /eV for the 2D band of unstrained 
graphene and the 2D3 and 2Di subbands of strained 
graphene at e = 2.6%, respectively. It has been argued 
that this strain might be about three times larger than 
the actual one^ As shown above, our strain rates over- 
estimate the experimental ones of Ref. [l^ by « 100%. 
Therefore, we compare the measured rates to the calcu- 
lated ones of 95 and 99 cm _1 /eV for strained graphene at 
e = 1.3%. Our rates for unstrained and strained graphene 
underestimate the experimental ones roughly by 5%. The 
clarification of the origin of this discrepancy requires fur- 
ther theoretical and experimental efforts. 

Within the simple model of linear electronic and 
phonon dispersions in the vicinity of the K point, the dis- 
persion rate has been expressed as {2/h){v p h/vF), where 
v p h and vf are the phonon group and electronic veloc- 
ities near the K point, respectively, and fi is Planck's 
constant over 2irJ£ In the case of tensile Z-strain, it has 
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FIG. 5. Scaled Raman shift of subbands 2Di and 2D 3 as 
a function of El at different e and at 9 — 0° and 30°. The 
triangles are the shifts of the 2D band of unstrained graphene. 
The lines are guides to the eye. 



been argued that vf decreases with increasing strain, 
the decrease being larger along path 1 than along paths 
2 and 3. Based on these arguments and the assump- 
tion for isotropic phonon group velocity, it has been con- 
cluded that the expression for the dispersion rate above 
describes the observed behavior of this quantity^ On 
the other hand, in Ref. [TH and HH, it has been demon- 
strated that Vph of strained graphene is anisotropic. The 
NTB calculations on Z-strained graphene show that, with 
increasing tensile strain, Vp remains almost unchanged 
and isotropic (the Dirac cones are displaced almost un- 
deformed), while v p h increases and becomes anisotropic. 
Moreover, the latter increases more for path 1 than for 
paths 2 and 3. The NTB velocities, substituted in the 
expression for the dispersion rate, yield within a few per- 
cent the theoretical one in Table III. The fact that con- 
tradicting theoretical arguments predict the same exper- 
imental observation indicates that more theoretical work 
is needed to reveal the correct dependence of the Fermi 
and phonon group velocities on strain. 



D. Dependence of the Raman shift of the 2D 
subbands on 9 



Figure 6 presents the Raman shift of the 2D Raman 
subband, extracted from the calculated spectra. For all 
sets of values of e, 9, and E L , the shift of the 2D2 subband 
is between that of the 2Di and 2D 3 subbands. At 9 = 0° 
the shifts of the 2D2 and 2D 3 subbands are equal and 
these subbands contribute to the composite band 2D2 + 
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e is represented by a 2 x 2 matrix with nonzero diago- 
nal elements e xx — e and e yy = ~ve (y is the Poisson's 
ratio). The tensor R is represented by a 2 x 2 matrix 
with elements ±sin(t9) and cos((9). The dynamical ma- 
trix for the strained structure can be expanded in a series 
in e, which consists of terms of even degree in sin(t9) and 
cos(t9). Therefore, the phonon frequency can be written 
as a polynomial of these trigonometric functions of even 
degree. In Ref. [3, a single function sin 4 (6 + 6>o) (#o is a 
proper constant phase shift, depending on the subband) 
was found sufficient to achieve a fair fit to the Raman 
shifts of the 2D subbands at e = 1%. For the considered 
sets of values of e and El , a single trigonometric function 
gives an unsatisfactory fit to the 2D subbands, while a 
fair fit requires a polynomial of sin 2 (6> + 6q) of at least 
fourth degree. Here, we did not attempt to fit such a 
polynomial to our data because the data points are not 
sufficient for the accurate determination of the polyno- 
mial coefficients. For the practical use of our results, an 
interpolation between the calculated data points should 
be done. 



FIG. 6. Scaled Raman shift of the 2D subbands (symbols) 
as a function of 8 at different values of e and El- The three 
branches at a given strain magnitude are the shifts of sub- 
bands 2Di, 2D2, and 2D3 in order of increasing (decreasing) 
shift at e < (e > 0). The horizontal lines are drawn at the 
Raman shift of unstrained graphene. 

2D 3 . Similarly, at 6 = 30° the shifts of the 2B 1 and 
2D2 subbands are equal and they should be observed as 
a single composite band 2Di + 2D2. At small El, the 
splitting of subbands 2D2 and 2D3 is small enough to 
make them difficult to separate in Raman experiments. 
We also notice that the curves of points for e = — 1% 
and e = 1% are almost centrosymmetric. This can be 
explained by the fact that a small tensile Z-strain de- 
forms the atomic structure of graphene similarly to a 
small compressive A-strain and vice versa. 

The derived shifts form specific patterns, which can be 
fitted with a polynomial in sin(t9) and cos(t9). The origin 
of such a dependence can be found in the fact that, upon 
a small uniaxial deformation, any vector x transforms 
into the vector (I + e')x, where e' is the strain tensor 
at strain angle 9. This tensor can be obtained from the 
strain tensor e at 6 = 0° by a rotation 9 by means of an 
orthogonal transformation, defined by an orthogonal ten- 
sor R: e = R~ 1 eR. 22 In the chosen coordinate system, 



IV. CONCLUSIONS 

We have presented a detailed study of the 2D Raman 
band of strained graphene within a non-orthogonal tight- 
binding model. We have obtained the dependence of the 
Raman shift of the 2D subbands on the strain magnitude 
and direction, and on the laser photon energy. The de- 
rived red (blue) shift of the subbands systematically over- 
estimates the available experimental data, which can be 
cured by using a constant scaling factor of 0.8. The shift 
rate and the dispersion rate of the 2D band change by 
up to 50% with varying the strain magnitude and direc- 
tion, and laser photon energy. This result underlines the 
importance of the theoretical and experimental investi- 
gation of the behavior of the 2D Raman band of strained 
graphene under different experimental conditions. Our 
results for the scaled Raman shift, shift rate, and disper- 
sion rate of the 2D band can be used together with the 
predictions for the shift rate of the G band 1 ^ for monitor- 
ing strain in graphene by means of Raman spectroscopy. 
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